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FREQUENCY STANDARD STABILITY FOR DOPPLER 
MEASUREMENTS ON-BOARD THE SHUHLE 

1.0 SUMMARY 


Attention is focused on crystal and atomic frequency standards, 
as elements of a system configuration for making one-Way Doppler mea- 
surements on-board the Shuttle. Emphasis is placed on crystal stan- 
dard performance in the Shuttle environment because crystal standard 
specifications have been used for the Shuttle frequency reference. 

g 

Sensitivities in the order of 1 part in 10 are noted for the 
effects of acceleration and vibration on crystal standard frequencies. 

9 

Typical aging rates from 1 part in 10 to values lower than 1 part in 
10^^ are considered to be feasible as a result of literature studies. 
Reports on the performance of spaceborne crystal standards show the 
effects on aging rates to vary. In one case, the rate was observed 
to be relatively unaffected. In another case, the rate was seen to 
vary from a positive value, through zero to a negative value, and 
finally to drift slowly back toward zero. 

Effects of the space environment on aging rate appear to include 

an effect of radiation. One report notes a fractional frequency shift 
-9 

of about 0.5 X 10 during a period of about one day. It appears that 
launching forces and radiation may combine to produce effects on the 
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frequency stability of a crystal standard which can be bounded in 
magnitude, but not precisely predicted at this time. 

Frequency stability terminology is reviewed in relation to the 
Shuttle application. Examples are provided from calibration measure- 
ments that are now being made, in preparation for tests on one-way 
Doppler equipment. These examples illustrate an approach for the 
development of an error model that would include stability effects 
of the ground station and Shuttle frequency standards. 
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2.0 INTRODUCTiai 


Tv/o frequency standards, or reference frequency sources, are to 
be used in the system configuration for making one way Doppler mea- 
surements on-board the Shuttle. The system error budget will include 
contributions from each of these sources. The reference frequency 
source that is to be used on-board the Shuttle vnll be subjected to 
rather severe environmental conditions. 

Consideration has already been given to the selection of a ref- 
erence source for the Shuttle baseline and a specification has been 
prepared^^^. The decision v/as made to use a crystal oscillator, 
rather than an atomic standard for this purpose. Atomic standards 
will be used, hov/ever, at the ground station installations. This 
report has been prepared to describe the short and .q term stability 
characteristics of crystal and atomic standards. Emphasis is placed 
on crystal oscillators because of the selection v/hich has been mads 
for the Shuttle baseline and the lomplexities v/hich are introduced 
by the Shuttle environment. 

Attention is given, first, to the definitions of stability and 
the application of these definitions to the Shuttle system and its 
mission. Data from time domain measurements are used to illustrate 
the definitions. These data are of interest because they are from 
initial calibration tests of a configuration that is to be used in 
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tests on one-way Doppler (OWD) equipment. Statistical treatments of 
these data illustrate methods to be used in error model development 

(2 3 ) 

and in the implementation of OWD test plans' * . 

The results of a literature survey to determine environmental 
effects on frequency reference sources is presented next. Finally, 
methods of standard frequency dissemination over radio frequency 
carriers are noted as a possible means of measuring absolute accu- 
racy and long-term stability characteristics during tests on OWD 
equipment. 

The scope of this report is limited to the results of a litera- 
ture survey and initial calibration tests on a configuration for eval- 
uation of OWD techniques and equipment. It should serve as back- 
ground information and one in a series of steps tov/ard the develop- 
ment of error models and performance predictions for one-v/ay Doppler 
measurements on-board the Shuttle. 
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3.0 STABILITY TERf-lIHOLOGY AND DEFINITIONS 


The stability of a frequency standard is usually described by a 
statement of long term stability and a statement of short term sta- 
bility. These descriptions suggest time domain definitions. Other 
descriptions refer to phase noise sidebands, suggesting frequency 
domain definitions. These tv/o domains are related, mathematically, 
and are equally useful in a general description of frequency stan- 
dards and their stability. In a particular application, hoirever, 
one definition may be more useful than another because it has a 
more direct relation to a measurement method or an observable 
characteristic. 

Several authors have provided helpful discussions^^”^^^ of 
stability definitions and terminology. Some of these discussions 
emphasize the importance of short-term stability because of their 
application to tivo-v/ay Doppler systems in v/hich a received signal 
can be compared, after a relatively short period of transmission, 
v.'ith the original source v/hich determined the frequency of the 
signal. In other applications, such as time-keeping, long-term 
stability is more important. Both short-term and long-term char- 
acteristics should be considered in an evaluation of standards for 
use in making one vray Doppler measurements on-board the Shuttle, 
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One of the most basic terms that is used in discussing the sta- 
bility of a frequency source is "frequency deviation" (Af). This is 
the difference between the source frequency at a particular instant 
(f|.)» or averaged over a knov/n interval of time, and the nominal (f) 
or mean frequency (f^jj) of the source, fleasured or computed frequency 
deviations are often normalized v/ith respect to the nominal or mean 
frequency and then referred to as fractional frequency deviations (F). 



Both long-term and short-term deviations can be illustrated by 
the use of Figure 1. Units of time and frequency are indicated in 
Figure 1, but are not defined. This permits the use of the illus- 
tration v/ithout regard for a strict line of demarcation betv'een 
long-term and short-term deviations. Consider, now, this illustra- 
tion to represent the frequency reference on-board a Shuttle. At 
launch, the frequency has the initial (f^) value indicated. The 
frequency has a relatively short-term fluctuation and an upv/ard 
drift v.'hich can be called a long-term instability. After t^ inter- 
vals of time, a long-term average freouency can be identified. If 
the long-term drift has been investigated for this particular ref- 
erence and has been found to be reasonably oredictable, this drift 
could be removed by comoutations within the navigation computer. 
Alternatively, periodic corrections could be made to the oscillator 
frequency throughout the flight. In either case, the deviation of 
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FRIiQUHNCY IN 
ARBITRARY UNITS 



TIMO IN ARBITRARY UNITS 





the trend line from the initial frequency line, as noted at t 2 » 
could be corrected by adding a predictable frequency change to the 
initial frequen^jy. 

Throughout the Shuttle flight, the actual frequency will fluc- 
tuate around the trend line. To discuss this fluctuation, it will 
be convenient to define a deviation (d 2 ) from the trend line. A 
statistical description of this deviation will be useful in pre- 
dicting Doppler measurement errors if a trend line computation or 
correction is made on-board the Shuttle. Alternatively, a long- 
term average frequency might be predicted for a flight. In this 
case, a statistical description of the deviation (d^) of the actual 
frequency from the long-term average frequency v/ould be needed. In 
the simplest case, the initial frequency would be used as the ref- 
erence throughout the mission. A statistical description of the 
deviation (d^) from this initial frequency would then be used to 
estimate Doppler measurement errors from this source. 

It will be convenient to have an expression for frequency as a 
function of time, f(t), for a particular source. Acceptance of a 
straight line approximation to the long-term drift permits the 
following expression to be v/ritten. 

f(t) = ^0 + + 1^ (2) 
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Th’ slope of the trend line in figure 1 is introduced as and n^ 
is the noise-like fluctuation of frequency, around the trend line. 

Tlie random nature of n^ prevents the writing of a completely deter- 
ministic description for the frequency as a function of time. The 
fluctuating frequency tem (n^) can, hov/ever, be described in terms 
of statistical parameters which can be estimated from laboratory 
measurements. 


3.1 A TIME DOMAIN DESCRIPTION OF SHORT-TERM 
FREQUENCY STABILITY. 


Some frequency standards have a negligible long term deviation 
(represented by d^ and d^ in Figure 1) and others have this quality 
over relatively short time intervals. Consider, then, either of 
these ti70 classes of oscillators and the need to define the Shuttle 
frequency reference source over a relatively short interval of time. 
During this interval, Doppler cycles are to be counted. For the 
time interval to be used here, the long term drift (dp dg) can be 
neglected. Alternatively, it can be assumed that this long term 
component of deviation is correctly oredicted and removed. The 
short term deviation (d 2 ), averaged over the period of interest is 





d2 dt 


(3) 
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It is assumed here, for convenience, that the counting interval is 
from t 2 to t^ (Figure 1). The rms value of the short term devia- 
tion, over this counting interval, is 



Equations 3 and 4 can be used to define average and rms fractional 

deviations (F„., F,,„, ) by dividing Af „ and Af by the nominal 

or mean frequency (f). 

Time domain descriptions are also v/ritten^^^ as a function of 
the phase noise tern in the following expression for the signal, 
e^Ct), from a frequency standard. 

SjCt) = [A + a(t)] cos [ojt + v>(t)] (5) 

The amplitude mean and amplitude fluctuations (amplitude noise) are 
represented by A and a(t), respectively. Amplitude fluctuations 
are not of first order significance in this consideration of freouency 
and phase stability and will not be given further attention at this 
time. 


The bracketed phase term in Equation 5 is a function of time (t) 
and the nominal frequency (f), since w = 2-nf. It is also a function 
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of the short-term frequency fluctuations, A theoretical ex- 

pression for instantaneous fractional frequency deviation. 


y(t) 


0 ) * 


( 6 ) 


has been used^^^ as a basic time domain description. A sampling pro 
cedure is generally used to obtain a practial measurement of y(t) 
and the sample collection requires an observation time (t). An aver 
age value for y(t), over the period of measurement is written as 

k+T 

y(t) dt (7) 

for kth period of x seconds duration. A formal definition for sta 
bility is the infinite time average of the sample variance of tv/o 
adjacent averages of y(t). This definition called the 




Allan variance, implies an infinitely long data record. However, 
a good estimate of stability is generally considereo to be obtainabl 
by use of a limited data record of m sample differences. 


J 


1 

- E 

m 


^^k+1 ■ ^k^^ 


k=l 


(9a) 
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This sum of squared deviations is obtained by subtracting the first 
satrr''e from the second, then the second from the third and continuing 
in this manner until m deviations are processed. Another expression 
(9b) for Oy (t) pairs the first and second samples, then the third 
and fourth, using each sample in only one deviation compu .tion. 

1=1 

A more general expression (9c) is sometimes used^^’®^ for (t). 
Whereas (9a) and (9b) are formed bv use of two adjacent samples, (9c) 
uses a sequence for which N is greater than 2. This is a more 
general form of the Allan variance. 

(N. T. T>= (7 - I ^ 7^)^ (90 

3.2 FREQUENCY DOMAIN DESCRIPTIONS AND THEIR 
RELATION TO THE TIME LOMAIM 

The work of standardization committees and references which 
(7 B) 

propose standards' ’ ' refer to the frequency domain description as the 
"first definition of frequency stabiliti'". The time domain descrip- 
tion is called "the second definition of frequency stability". The 
"’’rst definition of frequency stability (frequency domain) is said 
to be the one-sided spectral density of the fractional frequency 
fluctuations, expressed on a per hertz basis. Although it is 
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theoretically possible to do so, an ac' rate is not generally des- 
cribed in the frequency domain. Other types of noise, hov/ever, such 
as thermal noise, shct noise, and flicker noise, are readily described 
by frequency domain expressions for the one-sided, noise spectral 
density. 

( 8 ) 

A similarity has been noted' ' bet\/een some of the methods used 
to r asure stability in the frequency domain and methods used to make 
stability measurements in the time domain. A spectrum analyzer is 
used, of course, to provide a frequency v/indov for frequency domain 
measurements. A gated counter provides a time windov/ for time domain 
measurements. Generally speaking, stability measurements are made by 
use of phase or frequency discriminators, by use of a spectrum analy- 
zer or a ivave analyzer, or by use of a counter. Phase detector re- 
cordings and records of cycle counts ever knov/n time intervals pro- 
vioe the most direct measure of stability in the time domain. In 
either of these cases, the measurement may bo interpreted as a measure 
of the random phase noise, <r(t), in Equation 5. 

The autocorrelation function of the phase of a ‘<ave is defined 
to be an infinitely long time average, 

Ry,(T) - lim % j <fi{t + t) v(t) dt, (10) 

J 00 Jq 
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v/hich is adequately approxinated^^^ by a computation from a suffi- 
ciently large nunter of measurements. 

(t + x)^(t) (11) 

Since autocorrelation functions and power spectral densities are 
Fourier transform pairs, the phase spectral density can be written as 


Sy,(f) = /” M-f) e V ^ ^ 

wOO 


( 12 ) 


The corresponding one-sided spectral density is 

^00 

s^(f) = 2 I Mt) cos (z-rrf t) dx ( 13 ) 

•'o 

Differentiation in the time domain corresponds to multiplication 
by jZtrf^ in the frequency domain. An expression for frequency 
power density is obtained, then, as 

S^(f) = (2TTf^)2 S^(f) (14) 

This power spectral density of the frequency fluctuations is normalized 
to obtain the pov.'er spectral density of the fractional frequency fluc- 
tuation v/ith respect to the nominal frequency. 

= -V Mf) 

f^ 


F = ^ ^f) 


(15) 



Fourier transformations will, then, permit time domain measure- 
ments and relations to be described in the frequency domain. Conversely, 
frequency domain data and relations can be described in the time domain. 
In spite of the equivalence between time and frequency domains. Reference 
8 states "in order to have a comprehensive and sufficient measure of 
frequency stability, it is preferred that the measurements involve both 
frequency and tine domain techniques". Time domain averages of 1 to 100 
seconds correspond to 1 to 0.01 Hz offsets from the carrier in the fre- 
quency domain. Special attention must be given to spectrum analyzer 
resolution vrfien measurements are to be made in this region. Counters, 
however, are easily gated for this range of time averaging. Because time 
domain measurements v/ill be easier to make over this range of averagino 
periods, emphasis will be placed on this aoproach at this time. 

3.3 ^€ASUREMENTS WITH A TIflE DOMAIN 
CONFIGURATION 


The configuration of Figure 2 has been used to obtain stability 
data for several oscillators with different stability specifications. 

The vector voltmeter (HP 8405A/ measures the amplitudes and phase differ- 
ence for two input signals. An output voltage that is normally used for 
recording the phase difference is, in this application, used to start and 
stop the period counter (HP 5245L). If the two frequencies (f-| and f^) 
were exactly the same, the phase difference would not change with time. 
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Fl}>uro 2. - A configuration for measuring frequency 
difference periods and tlie period fluctuations. 



In this case, the voltage at the input to the period o^unter would be 
constant and the counter would not start or, if it were counting, it 
would not stop. The usefulness of this configuration depends, there- 
fore, on the existence of a small difference betv/een the two oscillator 
frequencies. 

Consider, as an example, that the oscillator frequencies differ by 
0.1 Hz. The period of the phase difference versus time wave is, in this 
case, 10 seconds. The counter clock frequency is, nominally, 1 MHz. 

The phase difference (0^) between the two oscillator signals increases 
at a constant rate, or decreases at a constant rate, if the two frequencies 
are constant. 


®d ~ 27r(f2 - f-|)t (16) 

The vector voltmeter is designed to provide an increasing voltage when 
the frequency of the input to Channel B is higher than the frequency of 
the input to Channel A. Full scale indication results for a phas? differ- 
ence of plus or minus 180 degrees. The voltage ramp increases from a 
minimiBn value for = -180° to a maximum value for 0{j = +180°. It 
then drops abruptly to the -180° indication and continues increasinq at 
the constant rate. 

When the 9(j(t) wave initially drops from +180° to -180°, the period 
counter is started. On the next drop, the counter is stopped. The 
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count that is accumulated Is 


N = 





(17) 


and is measure of the perioa, Tp. During this period, the phase differ- 
ence has changed by 360° or 2ir radians. From a rearrangement of (17), 
then, 

(f2-fi) = — (18) 

The measurement of N, as defined by (17), is, nominally, in microseconds. 
The accuracy of this measurement is dependent upon the accuracy of the 
clock that generates the f^ frequency. A sumary of the stability specifi- 
cations for the HP 5245L clock are included in Appendix A. As noted 
there, a measurement over a period of 10 seconds should not be affected 
by the aging rate, but short term frequency fluctuations v/ill normally 
occur within this period. A temperature change of one degree vrill cause 
a systematic frequency shift that v/ill be comparable to the short term 
frequency fluctuations. Line voltage changes of several percent would 
cause systematic changes of this same order of magnitude. Collectively, 

9 

these effects could cause a frequency change of about 1 part in 10 
v/hich is 0.001 Hz for the 1 fiHz clock. Tlie net effect for a 10 second 
measurement could be 0.01 Hz. Since this is beyond the resolution of the 
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8 digit display, it will not affect the measurement. Similar reason- 
ing leads to the conclusion that a clock frequency inaccuracy of less 
than 0.1 Hz will not affect the measurement for a 10 second oeriod. 

The normal error of one count will, however, cause an error in the 
eighth digit. 

A typical period measurement of = 8.115878 seconds was made 
vnth the configuration of Figure 4, using a synthesizer driver (HP 511 OB) 
as a source of f-| and the counter clock as a source for ^ 2 - 
counter clock also served, of course, as a source for f^. The fre- 
quency difference between the synthesizer driver and the counter clock, 
averaged over the 8+ second period v/as, then, 0.1232153 Hz. The next 
measurement v/as 8.119960 seconds, giving an average frequency difference 
of 0.1231533 Hz. The difference in the t\r/o frequency measurements is 
0.000062 Hz. This is a 6 part in 1 0 ^ Vl uctuati on in the frequency 
difference betv/een the two oscillators. The measurement does not 
assign the fluctuation to either f-j or f^, but to the combined 
fluctuation of f-j and f^. As noted in Appendix A, this fluctua- 
tion is within the specification of the counter clock, but exceeds 
the specification of the synthesizer driver. As a first approximation, 
the fluctuation could be attributed to tho oscillator that is knov;n to 
be less stable. Or, it the tv;o oscillators are approximately equal in 
stability, the fluctuation could be assigned to both of them, on an 
equal basis. A better evaluation can be had, hov;ever, by the collec- 
tion of more data and by the use of statistical methods of analysis. 
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4.0 TIME DOMAIN DATA ANALYSIS AND EVAL-, \TIOfl 


The Doppler measurement is to be made as a count of cycles that 
represent the one way Doppler shift. The cycle count will be made 
over some finite period of time that v;ill be used v/ith the cycle count 
to estimate an average velocity or a range increment. By its nature, 
then, the Doppler measurement is a r«asure of frequency, averaged over 
a period of time. Consecutive . counts of Doppler cycles will be made 
and summed to provide an integral of the velocity, which is a sum of 
range increments. 

The nature of the inplenentation and application of Doppler mea- 
surements in the Shuttle operation involves a series of data points, 
arrayed in time. It is normal, then, to look first at a set of data 
presented as a time series. Table I contains ti.'o sets (x^ , X 2 , x^,..., 
x^g) of period measurements, made with the test configuration of Figure 
2. As noted to the right of each data point, the reciprocal of the 
measurement v/as taken to obtain the frequency difference between the 
fc/o oscillators that v;ere used in the test. The E+ symbol indicates 
the next operation, which is the entry of the reciprosal into a sub- 
routine that connuted the mean freouency diffe»e«ce (x) and the stan- 
dard deviation (s) for the set of 16 measurements. Next, the difference 
was obtained betv/een each measurement reciprocal and the mean frequency 
difference. Each of these deviations from the mean was divided by the 
nominal frequency of the t./o oscillators (10 Hz), to obtain a value of 
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TABLE I - PERIOD MEASUREMENTS MADE BY USE OF 
A CONFIGURATION LIKE THAT OF FIGURE 2 
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fractional frequency deviation. These fractional frequency deviations 
about the mean are plotted in Figure 3 for the data set on the left of 
Table I. The plot is linear in time and covers a total time of about 
4.3 minutes, a period of time that could represent an orbital pass of 
the Shuttle at some typical ground-track distance from the transmitting 
ground station. 

4.1 A PROBABILISTIC MODEL OF THE DATA SET 

The most basic features of the configuration of Figure 2 are iden- 
tical to those of the configuration that v.'ill exist during Shuttle 
operations. This is particularly true with regard to the test oscilla- 
tors of Figure 2, which can represent the ground station standard and 
the Shuttle reference. Simplification of the configuration, from that 
of the Shuttle and a ground station to the test configuration, places 
emphasis on the frequency sources and the errors that they will intro- 
duce in the Doppler measurements. Information from tests with this 
type of test configuration will be useful in predicting the accuracy 
of Doppler measurements on the Shuttle. If, for example, one of the 
test oscillators is a ground station frequency standard and the other 
is a Shuttle frequency reference, the configuration can be used to ob- 
tain an estimate of the error budget contribution from these two 
sources. 
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^-igurc 3. - A time- series jilot of 
J c V i a t i o n s f o r t wo c r v s t a I 


I iMct lonal freqiioncv 
'se 1 1 1 a t ors . 


A probabilistic model is needed for the fractional frequency fluc- 
tuations that are measured for the two oscillators of Figure 2 if these 
data are to be used in making accuracy predictions. Such a model is 
also important in data evaluation and in the development of a meaning- 
ful specification of oscillator characteristics. To obtain such a 
model, measured data were compared to a normal probability density 
function. The normal, or Gaussian, density function (y) was selected 
for first consideration because many naturally occurring noise phenomena 
have the characteristics of this function. The mathematical expression 
for this function. Equation 19, 

_ (x - u)^ 

y = — ==- e 2a (19) 

a\2iT 

i/as changed by use of the standard deviate (z) as the independent 
variable. Also, a was deleted from the denominator because dx = adz. 


- X - u 


( 20 ) 


The model used, then, on a trial basis, v/as 


y = 


r 


■^2ir 


( 21 ) 


Day-to-day variations in the mean (u) and the variance (a), and 

2 

estimates of the mean (x) and variance (s ), were removed from the 
model in this v/ay. 
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Nine sets of data, such as that on the left of Table I, i/ere used 
to compare the measured frequency fluctuations v/ith a normally distri- 
buted population. The comparison is illustrated in Figure 4. The 
computed mean (x) and standard deviation (s) for each set were used to 
convert the measurements (xp X 2 » x^.-.-.x^g) of the set to d corres- 
ponding set of standard deviates (Zp Z 2 , Z 2 ,...,Z^g) by use of Equation 
22 . 


7 


Each of the 144 values of Z obtained in this v/ay v/as assigned to a 
group that corresponds to an interval along the Z axis. The number 
of Z values* Mn each interval i. represented in Figure 4 l»y a vertical 
bar. The bar v/idth represents the Z-axis interval. The smooth curve 
of Figure 4 is a plot of Equation 21, normalized to peak v/ith the 
experimental frequency distribution. 

The comparision provided by Figure 4 indicates reasonable agree- 
ment betv/een measured data and a normally distributed population. The 
most significant discrepancy appears to be around z = -0.8. The use 
of more data might re 'jce the differences betv/een the experimental 
and theoretical distributions. On the other hand, the experimental 
function may prove to be unsymmetrical. For example, some real physi- 
cal phenomenon might skew the experimental function toward the direction 
of the long-term difference in the drift rates of the two oscillators. 


4-6 




DIFFERENCE FROM THE MEAN IN STANDARD DEVIATE UNITS 

.re 4. “A comparison of measured data with a normal 

probability density function. 


Since the agreement between the theoretical model a. ' the experi- 
mental data is easonably good, a normal density function is accented 
as a reasonably good model for the short-term frequency stability. The 
theoretical plot of Figure 4 can be used to predict the probability 
of experiencing fractional frequency fluctuations of a particular amp- 
litude, For the data plotted, one standar*^ deviate corresponds to a 

_5 

standard deviation of 2.889 x 10 i!z, or a fractional frequency fluc- 
tuation 0* 2,889 X 10'^\ These, then, are one-sigma values. We pre- 
dict, therefore, that 68.26 percent of all future measurements of 
short-term stability, for these tv/o oscillators and these conditions, 
will not exceed a fractional frequency fluctuation of 2.889 x 10'^ \ 

Correspondingly, we predict that 99.73 percent of such measurements 

-11 

v/ill not exceed the three-sigma value of 8.666 x 10 
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4.2 VARIATIONS IN VARIANCE 


Three ways of computing the variance of a set of period, phase, 
or frequency measurements were noted in Section 3.1 and each of these 
may be preferrable in a particular situation. It is not clear, yet, 
which of these is best suited for the evaluation of Doppler measure 
merts in the Shuttle application. 

The standard deviation that is noted in Table I was computed in 
the conventional manner, as the square root of the variance, which is 
the sum of squared deviations from the mean, divided by the degrees of 
freedom. 


s = 


Z(x - x)^ 
n - 1 




1/2 


(23) 


The variance for the set of measurements on the right of Table I would 
be s^ = (1.316 X 10"^)^ = 1.733 x 10“^®. Using Equations 9a and 9b, 

the variance values for this same set of data are computed to be 

c/(t) « lij E -7,)^ - 1.549 X 10-’» 

and 

N 

H = 1.328 x10-'» 
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The standard deviations for these three computations are, respectively. 


s = 1.316x10'®, 
s = 1.244 X 10'®, and 
s = 1.152 X 10'®. 

It can be seen, intuitively, that long-term drifts will have less 
effect on the second and third calculations than on the first. In 
fact, the second and third calculations effectively remove long term 
drifts from the computed stability value, unless the yj^ measurements 
are rather long term averages. This may be desirable, particularly 
if a specifically defined short-term stability value is needed in 
the application under consideration. Long-term drifts are important, 
however, in the Shuttle application. In this application, any period 
greater than the sample period, and particularly the interval between 
sampling instants, must be examined for long-term stability character- 
istics. 


4.3 REGRESSIONS IN TIME 

Data from Table I can be used to estimate both short-term and 
long-term stability. A regression, in time, can be made within the 
set to obtain "long-term" stability characteristics for periods in 
ohe order of several times the sample duration, in this case about 
8 seconds. A regression, in time, using the means from several data 
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sets can be used for longer-ten stability estimates. The set spacing, 
in time, controls the period over which the stability is to be 
specified. 

To obtain an estimate of the per-day difference in drift rates 
for two oscillators, data were taken over a one month period, using 
the circuit of Figure 2. A set of 16 or more measurements, such as 
those in Table I, was made at approximately the same time on 16 diffe- 
rent days during the month. The mean (x) from each of the 16 sets of 
data was used in a linear regression, with time in days (t^) as the 
independent variable. In this way, an expression for the frequency 
difference (f^) between these two oscillators was computed to be 

fj, = 0.0605 + 0.00221 t^j Hertz. 

The stabilities for these two oscillators, as specified by the manu- 
facturer are 

lOlA: 7 X 10*^ Hz per day 
51108: 3 X 10~^ Hz per day 

Assuming these rates to be positive, as is most common for crystal 
oscillators, the difference in the drift rates is expected to be 

_3 

4 X 10 Hz per day. The measured drift difference is the computed 

_3 

coefficient of t^ in the regression equation, 2.21 x 10 Hz per 
day. 
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A correlation coefficient (r) can be computed as a part of a 
regression analysis, to indicate the degree of the relation between 
the two variables. For the test just described, the correlation 
coefficient was computed to be 0.4133. If the two variables were 
completely unrelated, or uncorrelated, the correlation coefficient 
would have been zero. If the data points had all been on the re- 
gression line, the coefficient of correlation would have been unity. 
The correlation coefficient is considerably less than unity. This 
indicates a need for a nonlinear regression model, a random scatter 
of the data around the linear model, or both of these conditions. 

The square of the correlation coefficient is sometimes called 
the coefficient of determination and is a mea jre of the extent to 
which the regression has explained the variance in the original se- 
quence of data-set means. The unexplained, or residual, variation 
is (1 - r ) times the original sum of squared deviations of x^ values 
from their mean. The original variance for the set of means was 
1.822 X 10’®. Of this amount, (0.4133)^ (1.822 x 10’^) = 3.11 x 10“^ 
was explained by the regression and is considered to be a 24 hour 

drift rate. The residual amount, [1-(0.4133)^] (1.822 x 10"^) = 

-5 

1.511 X 10 is still unexplained. As noted earlier, some of this 

2 

residual variance (s^ ) might be removed by a nonlinear model. Other- 
wise Doppler measurements made with these two oscillators would be 
in error by an amount that is proportional to the square root of this 
residual variance, (1.511 x 10“^)^^^ = 3.89 x 10"^. 
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Another useful application of this type of regression analysis is in 
the evaluation of oscillator warmup (recycle) time. A linear 
regression on a set of measurements for a 5 MHz oscillator, within 
an hour after it was turned on, resulted in a correlation coefficient 
of 0.9999, indicating a highly time-correlated frequency. On the 
following day, the drift rate had decreased by more than two orders 
of magnitude and the correlation coefficient had dropped to 0.3368. 
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5.0 EHVIRONf!EflTAL EFFECTS OF SHUTTLE FLIGHTS 


Environmental stimuli v/hich can be expected during Shuttle flights 
inci ude: 

(a) vibration, 

(b) shock, 

(c) acceleration, 

(d) temperature, 

(e) radiation, 

(f) supply voltage changes, 

(g) magnetic fields, and 

(h) electrostatic fields. 

Some degree of isolation can, of course, be specified for each of these 
stimuli. Such isolation adds cost and weight to the Shuttle design, 
however, and should not be over specified. It is important, then, to 
define the stimuli levels and compare these vn'th predictions of their 
effect on each type of frequency standard. 

Current plans for use of OWD measurenents on the Shuttle only re- 
quire such measurements f’urinn orbital flight phases. It is assumed 
here that these plans include the use of OHD data during orbital maneu- 
vers, such as orbit changes or corrections and rendezvous operations. 
Environmental conditions, then, are expected to involve relatively 
little vibration and shock during use of the OWD equipment. The equip- 
ment would, hov/ever, be required to '.n'thstand the conditions of launch. 
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re-entry, and landing without damage and without any effects which 
would degrade future operation. Temperature and radiation condi- 
tions are to be controlled, to some extent, during orbital flight 
and during the launch and landing phases. Levels of the other 
stimuli noted above will probably not depend greatly on the flight 
phase. 

Improvements in the navigation solution might be realized during 
launch and return to earth, as v/ell as during ort)ital flight by Doppler 
measurements. For this reason, the effects of greater shock and vibra- 
tion levels during launch and the re-entry through roll-out operations 
are of considerable interest. This interest, along v/ith the possibility 
that high levels of shock and vibration may cause a permanent degrada- 
tion of data or an effect that has a long recovery time, focuses signi- 
ficant emphasis on launch and re-entry conditions. 


5.1 PERFORf^ANCE AND TEST REQUIREMENTS 

Accuracy and stability specifications on the Shuttle frequency 

reference are significantly influenced b'/ the system requirements for 

(14) 

one-v;ay Doppler measurements. Table IP identifies these system 
requirements, v/hich include effects of ground station frequency accu- 
racy and stability, propagation, receiver noise. Shuttle frequency 
accuracy and stability, and environmental conditions. 
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TABLE II - ONE-WAY DOPPLER SYSTEM ACCURACIES 


Counting Time 
(Seconds) 

Maximum Noise Error 
(Cycles, 3a) 

0.1 

0.3 

0.5 

0.3 

1.0 

0.4 

2.0 

0.5 

10.0 

2.3 

30.0 

6.0 

I 

60.0 

12.0 

600.0 

120.0 


NOTE: Maximum Allowable Bias Error is 2 Hz for a 1 day 
period, at the S-Band frequency. 
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The effects of environmental conditions on the Shuttle frequency 
source are estimable through the use of stability factors that relate 
frequency stability to each environmental condition, along with a de- 
finition of the environmental conditions. Table summarizes 

the environmental conditions which have been specified as requirements 
for the Shuttle frequency source. Figure 5 defines the random vibra- 
tion specification for the Shuttle as a function of frequency. Table 
is a list of performance characteristics which must be met by 
the Shuttle frequency source. The Shuttle frequency source must meet 
the performance requirements of Table IV during exposure to, and after 
exposure to, any feasible combination of the environmental conditions 
noted in Table 

Compliance with the long-term stability requirement, noted in 
Table IV is to be verified^^^ by at least 100 frequency measurements. 
These 100 measurements are to be equally distributed over a 24 hour 
period. Each of the measurements is to be taken over a 100 second 
averaging time. 

The short-term stability of the Shuttle frequency source is to 
be verified^^^ for each of tne averaging times (0.5, 1.0, 2.0, 10.0 
seconds) noted in Table IV. At least 100 frequency measurements are 
to be made at each of these averaging times. The measurements for a 
particular averaging time are to be made at intervals that do not 
exceed one minute. 
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TABLE III - environmental CONDITION SUMflARY 


(References 1, 16 and 17) 


Environmental 

Condition 

Requirement 
Speci f i cati on 

Vibration 

11.78 g rms random vibration 
(see spectrum in Figure 5) 
for 2 hou."s along each axis. 

Shock 

Landinq, one application of a 
1.5 g (peak) pulse of 260 milli- 
seconds duration. Also, more 
applications of lower peak and 
longer duration. 

Acceleration 

5 g for 5 minutes in both 
directions along each axis. 

Temperature 

35°F to 130°F 

Radiation 

None 

Supply Voltage 

24 to 32 volts 

i 
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TABLE IV - PERFORMANCE REQUIREMENTS FOR 
THE SHUTTLE FREQUENCY SOURCE 
(Reference 17) 


I erformance 
Characteristic 

Requi rement 
Specification 

Frequency 

4.608 MHz 

Adjustment 

part in 10^ 

Resolution 

1 part in 10^^ 

Stability 


Lonij-Term 

9 

1 part in 10 or less over a period of 
24 hours 

Short-Term 

1 part in 10^^ (rms) or less for 
averaging times of 0.5, 1.0, 2.0 and 
10.0 s. 

Stabilization Time at 

Frequency error ar; ’ drift rate shall 

a temperature of 30° F. 
(after 72 hours off) 

not exceed 

1 hour 

+ 2 X 10"^ 

8 hours 

+ 7 X 10'^ 

24 hours 

+ 2 X 10’^ 

! 

Repeatability 

Within + 1 X 10“^ after 8 hours 

(after 48 hours off) 


Acceleration 

_o 

1 X 10 '7g or less 

Sensitivity 
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5.2 ENVIRONMENTAL EFFECTS ON CRYSTAL AGING 


One of the disadvantages of the crystal oscillator, insofar as 
the Shuttle applicction is concerned, is a long-term drift character- 
istic. This drift is generally called, or attributed to, aging. Dis- 
tinction is sometimes made of a related characteristic, called fre- 
quency stabilization^^®^ This characteristic, stabilization, mani- 
fests itself as a relatively rapid rate of frequency change, or drift, 
during the first few weeks or months after production processes are 
completed. 

The transition from a period of stabilization to the period 
characterized by the rate of drift called aging is gradual, rather 

M g\ 

than distinct. In fact, the drift rate history has been described' 
as an exponential lunction that can be extrapolated. The curves of 
'figures 6 and 7 illustrate this concept, following the shape and 
magnitude reported in Reference 19. The oscillators tested in that 

Q 

' tudy showed a drift rate of one part in 10 per month, within the 
first few months of testing. Data indicated, however, that a long 

g 

term rate of one part in 10 per year could be expected. 

Several of the environmental effects, to which space borne 

crystal oscillators are generally subjected, might causa frequency 

(20 21 1 

shifts or changes in the aging rate. Two observers' ’ ' have re- 

porte( aging rates for an oscillator on the Timation II satellite. 
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Figure 6. - Typica! aging behavior for a good crystal oscillator 
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Figure 7. - A linearized approximation to the lo 
drift of a good crystal oscillator' 



Prior to launch, the drift rate was about 1 part in 10^^ per day. 

During the initial period of operation in space, the rate went to a 
high of about 3 parts in 10^^ per day and followed an aging rate 
curve with a negative slope for more than 100 days. As noted in 
Figure 8, the aging rate oassed through zero and went as high as 
about 5 parts in 10^^ per day, in a negative direction. Then, it 
began to decrease again at a decreasing rate. The last data re- 
ported in the reference indicate a negative drift rate of about 2 
parts in 10^^ per day. A significantly different aging character- 
istic was reported fr(»n Timation I. Figure 8 shows a very low ard 
almost flat aging rate history for this oscillator. 

5.3 RADIATION EFFECTS ON CRYSTAL OSCILLATOR FREQUENCY 

(221 

Se~ (conductors are known ^ ’ to undergo changes when subjected 

to radiation of the type and levels experienced during space flight. 
Corresponding information regarding quartz crystals has also been 
found in the literature reviewed. A shift in frequency has been 
reported jjowever, to have been caused by transient radiation. 

The plot of Figure 9 shows a relatively abrupt change of about 4 parts 
in 10^^ in a negative direction. This change appears to have occurred 
over a period of about one day. The aging rate prior to, and after, 
the radiation occurrence appears to be in the order of 1 part in 13^^ 
or less. The change, then, was by a factor of about 40. 


5-n 



(U 

H 

2 

O 

z 

l-H 

CD 

< 





I L 

rr ro 



5-1 


DAYS Al-Tl-R LAUNCH 



o: 

o 

H 


U 

C/3 


Z 

o 

E-* 

< 


O 

iH 

O 




5-13 


Figure 9. ^ Frequency shift in a spaceborne crystal oscillator, 

attributed to rad iat i on (20) . 



According to 


( 23 ) 

' i the characteristic eftcc^* of radiation on quartz 
oscillators is the introd'ction of a negative drift component. This 
effect has been verified in the laboratory, and should be anticipated. 
The reason for the better performance of the Timation I oscillator, 
compared to the Timation II oscillator (see Figure 8), is believed to 
be the counteracting effect of radiation. A normally positive aging 
was almost nulled by the negative aging rate caused by radiation. Be- 
cause of this experience with crystal oscillators, two rubidium oscilla- 
tors art: to be flown on Timation III, along with one crystal oscillator. 
Radiation sensors are planned to obtain an estimate of the correlation 
between incident radiation and the aging rate. 

5.4 THE IMPORTANCE OF CRYSTAL MOUNTING DESIGNS 


In an effort to improve the frequency stability of crystal oscilla- 
tors, Warner^ investigated the effect of different mechanical mount- 
ing methods on the performance of an oscillator under conditions of 
changing temperature, shock, and vibration, .'he quartz crystal was 
plated with gold, mounted on its quiescent edge in an evacuated en- 
closure, and operated in a thickness shear mode. In this way, the 
mechanically vibrating portion of the crystal was isolated from the 
mounting. 

An initial step in the mounting design is the enclosure of the 
mount in an evacuated housing. The crystals developed and tested in 
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the Warner^ ' study were housed in an evacuated glass envelope, gen- 
erally cylindrical in shape. The maximum dimension of the envelope, 
the height, was 2 inches. The diameter of the cylinder was just under 
2 inches. Findings indicated an increase of 20 or irore in the Q if 
the crystal is housed by use of a vacuum-enclosed mounting. A max- 
imum Q value of 10 million was measured for shaped, vacuum-mounted 
crystals. The maximum Q of a flat quartz plate in air was about 
400,000. Plano-convex shaping was used to improve tanperature-coefficient 
control, with no apparent loss in Q. 

Three ribbon-shaped mounting elements rigidly supported the crys- 
tal with a minimum of thermally induced strain. This three-point 
mounting permits radial expansion of the crystal plate, but restrains 
the plate from translation or rotation during mechanical shock. With 
this new mount, the time for the frequency to recover to within a few 

p 

parts in 10 after a step in temperature was only 2 hours. The use of 
rod mounting, instead of the ribbon mounting, resulted in a 12 hour 
recovery time. Similarly, the time required for frequency stabiliza- 

g 

tion to about 1 part in 10 per day was about two weeks for the rod 
mounting and only two hours for the ribbon mounting. 

Static loads as high as 200 g's, in any direction, cause no 
apparent movement of the crystal with respect to its mounting plat- 
forms. Severe shock, such as a four-inch drop, will dislodge the 
mounting, however. No permanent frequency change greater than one 
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part in 10 was observed when 5 Mlz oscillators, using crystals with 
this mounting were subjected to 10 g vibration of 2000 Hz. 

Improved performance of oscillators in which this crystal and 
mount design was used is attributed to the relief of mechanical 
strain, including thermally induced strains. A gravity-induced strain 
was observed, however, to cause fractional frequency deviations as high 

9 

as 1.5 parts in 10 . 

5.5 VIBRATION, SHOCK, AND ACCELER-MION 

The effects of vibration, during the launch phase, and the shock 

and deceleration jduring re-entry, are expected to require attention 

in the design, specification, and testing phases of the Shuttle-borne 

frequency reference development. Crystal oscillators are reputed^^^^ 

-9 

to have a vibration sensitivity of about 10 /g. Such a sensitivity 
is generally attributed to oscillators that use ribbon mounted crystals 
wher more specific information is not available from test data. Parti- 
cular attention to the mounting technique has reduced the vibration 
sensitivity to about 10"^ ^/g. 

When the crystal oscillator is subjected to sinusoidal vibration, 
sidebands are produced just as though a sinusoidal modulating source 
had been introduced. The modulation index increases in a linear manner 
with vibration amplitude over the limits of vibration reported in 
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Reference 24. The modulation index decreases with frequency, however, 
indicating a response to the modulating source that is typical of a 
frequency modulation process rather than a phase modulation process. 

The results of tests reported in Reference 24 are illustrated in 
Figures 10 and 11. The crystal mount resonance frequency is a domi- 
nant feature of Figure 11. This resonance effect was unacceptable in 
the pulse-Doppler radar application for which the tests were performed. 
To overcome the problem, a flexible cable support was designed for the 
entire crystal oscillator. The resonant frequency of this support 
cable was about 10 Hz. Above this frequency, crystal vibrations were 
attenuated by the supporting cable to a level that decreased at a high 
rate with increasing frequency. The approximate transmissibility of 
the cable support is shown, as a function of frequency, in Figure 12. 

Crystal oscillators from three manufacturers were subjected to 
acceleration and vibration tests at the Johnson Space Center recently 
The oscillators were not designed for the environmental condi- 
tions that have been specified for the Shuttle, but the performance 
measured is significant because it identifies the improvements that 
must be provided by more rugged designs. 

Static performance measurements were made before the oscillators 
were subjected to the vibration and acceleration tests. Data from 
these measurements are summarized in Table V. Tables VI and VII 
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Figure 10. Phase deviation of a crystal oscillat 
funct ion of sinusoidal vibration amplitude 





frequency for sinusoidal vibration 








TABLE V - TEST OSCILLATOR PERFORMANCE 
PRIOR TO ENVIRONMENTAL TESTING 
(Data from Reference 17) 
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summarize the acceleration and vibration test results, respectively. 

As noted in Table IV, two of the oscillators met the Shuttle specifi- 
cations for aging rate (1 part in 10^, or less, per dav) and short 
time stability (1 part in 10^^, or less, averaged over a one second 
period). The third oscillator met the short time stability, but 
failed to meet the aging »'ate requirement. 

Only one of the oscillators met the requirement for sensitivity 
to acceleration (1 part in 10^/g, or less) in both directions along 
each axis. A second oscillator was well within the sensitivity speci- 
fication on two of the three axes. The plots in Reference 17, of 
frequency shift as a function of acceleration, are nearly linear. 
Because of this, the possibility of modeling the effect and using a 
correction subroutine in the navigation program was suggested. 

Vibration test results show a need for a vibration isolation de- 
sign. None of the oscillators stayed within the required one part in 

Q 

10 during the vibration test. One of the oscillators did, however, 
stay essentially within the short-term stability requirement of one 
part in 10^® throughout the vibration test. It appears that the 
vibration amplitude would have to be attenuated by about 26 dB to 

g 

ensure a shift of no more than one part in 10 during vibration. If 
operation is not required during vibiation, isolation requirements 
would be greatly reduced. In this case, only the offset after vibra- 
tion would be considered in a vibration isolation specification. 
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TABLE VI - SUMMARY OF ACCELERATION TEST RESULTS 
(Data from Reference 17) 


Axis Along Which 
The Force of 
Acceleration Was 
Di rected 

Sensitivity, Parts in 10^/g 

Austron 

1120-1 

Freq. Elect. 
2037B 

Green ray 
YH 855-1 

+ X 

-0.1 

+1.0 

+4.0 

- X 

+0.3 

+0.8 

+2.4 

+ y 

+2.1 

-0.7 

-0.4 

- y 

f2.3 

1 

O 

-1.1 

+ z 

-0.5 

+0.8 

-1.3 

- z 

-O.'' 

+0.6 

-3.3 
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TABLE VII - SUIWARY OF VIBRATION TEST RESULTS 
(Data from Reference 17) 


Axis of 

Stability, Parts in 

10® 

V ibration 

Austron 

1120-1 

Freq. Elect. 
2037B 

Greenray 
YH 855-1 

Excursion 

peak-to-peak 




X 

20 

15 

★ 

y 

7 

15 

★ 

Z 

8 

7 

* 

Offset After 
Vibration 




X 

9.4 

1.0 

* 

y 

1.0 

1.0 

* 

z 

i.O 

0.0 

* 

i Short Term 
1 peak 




X 

0.6 

0.11 

* 

y 

0.3 

0.09 

* 

z 

0.3 

0.05 

if 


* This oscillator did not survive the vibration test. 
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5.6 ENVIRONMENTAL EFFECTS ON ATOMIC 
FREQUENCY STANDARDS 

The most significant advantage of atomic frequency standards 
over a crystal standard, in the Shuttle application, is the improved 
'.ong-term stability. The sunmary of frequency standard characteris- 
tics in Appendix 6 shows this advantage to be more than one order of 
magnitude for a rubidium standard. An even greater advantage is noted 
there for a cesium standard, for which the long-term drift is essen- 
tially undetectable. Little consideration need be given to other 
atomic standards because of size and weight disadvantages, or because 
of current developmental status. 

Size, weight, and complexity comparisons of Appendix B are rather 
decisive factors, causing the cesium standard to be deleted from com- 
petition for use on-board the Shuttle. The cesium standard is of 
interest, however, in the evaluation of overall accuracy, which must 
include the effects of ground station frequency standard stability. 

The rubidium standard is seeii to introduce additional environmental 
effects, along with increased complexity. The increased complexity 
appears to cause little change in size and weight, compared to a 
crystal oscillator. It would also jffect reliability, however, to 
soma extent. Environmental effects of ambient pressure and magnetic 
fields are introduced, but are seen to be rather small with respect 
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to the long-term drift specification of one part in 10®. If future 
consideration is given to the use of a rubidium standard on-board the 
Shuttle* these sensitivities and the related ambient specifications 
should be used to predict error budget components. 

In general, environmental sensitivities which have been observed 
in both crystal and rubidium standards are lower in rubidium standards 
than in crystal standards. Sensitivitv to acceleration, tor example, 

Q 1 A 

is about 1 part in 10 (or perhaps 1 part in 10 for special designs) 
for the crystal standard and less than 1 part in 10^^ for the rubidium 
standard. This implies a lower sensitivity to vibration, but data 
have not been found to support this conclusion. Temperature sensitivi- 
ties appear to be in the same order of magnitude, with a possibility of 
some advantage in the direction of the crystal oscillator. Radiation 
effects of about 2 parts in 10^^ for rubidium and about 4 parts in 10^^ 
for crystals have been reported by the same experimenters, but under 
considerab''y different circumstances. Environmental effects oii rubi- 
dium standards are described more completely in Appendix B. Also noted 
there is a relatively low warmup time for the rubidium standard, a 
possible advantage in the Shuttle application. 
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6.0 TIME AND FREQUENCY TRANSFERS BY RF LINK 


Discussions in previous sections of this report describe feasible 
methods of investigating short-term stability and note the importance 
ef both long-term and short-term stability. The absolute accuracy of 
the Shuttle reference frequency is also of importance, since an offset 
of this reference source, with respect to a ground station, will appear 
as a bias error in the Doppler measurement. The frequency of each 
ground station must, therefore, be accurately known. 

The transfer of an accurate frequency reference by way of a radio 

link provides a method of comparing or controlling a secondary standard, 

(251 

or reference frequency at a remote location. Table VHP summarizes 
the development of techniques for transmitting standard frequencies 
over radio links. Early broadcasts were made in the high-frequency (HF) 
band (3 to 30 MHz). Later, broadcasts were begun in the low frequency 
(LF) and very low frequency (VLF) bands because propagation conditions 
are more stable. In addition to the NBS media and -nethods noted 
in Table VIII, Loran C navigation transmitters and er VLF stations 
that are operated by the Navy provide a reference tc the Naval Obser- 
vatory time and frequency standards. Dissemination by satellite has 
also proven to be useful in recent tests. 

^251 

The NBS standard uses a long-beam cesium resonance tube' . 

12 

This standard has an estimated 3-sigma accuracy of 5 parts in 10 . 
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A working ensemble of NBS clocks are referred to this standard and, 

in turn, provide frequency references and control for NBS radio 

broadcasts. The WWV broadcasts are controlled to a long-term aver- 

12 

age frequency that is accurate to within 5 parts in 10 . Time sig- 

nals, as broadcast by WWV are accurate to within 5 microseconds. The 
NbS station in Hawaii, WWVH, has corresponding accuracies of 5 parts 
in 10^^ and 20 microseconds. VLF (3 to 30 kHz) and LF (30 to 300 

kHz) stations that are operated by the NBS, WWVL (20 kHz) and WWVB 

12 

(60 ’^Hz), are accurate to within 6 parts in 10 . This maximum 

degradation is said to occur during periods of phase correction. 
Broadcast accuracies, as noted here, denote deviations with respect 
to the standard NBS III. 
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TABLE VIII - SIGNIFICANT EVENTS IN THE 


DEVELOPMENT OF RF TIME AND FREQUENCY TRANSFER MEDIA 
(References 25, 26 ^nd 27) 


1 

Date 

Event 

1923 

UWV established to disseminate standard 
frequencies 

1935 

Time pulses added to WWV broadcasts 

1945 

Tirae-of-day voice announcements added to WWV 
broadcasts 

1948 

WWVH placed in operation to improve coverage 

1956 

WWVB began operation on 60 kHz 

1960 

WWVL began operation on 20 kHz 

Jan 1960 

All NBS broadcasts have been referenced to a 
cesium standard since this date 

August 1962 

First satellite time experiments with Telstar. 
(Between USNO and Royal Greenwich Observatory 
in England) 

February 1965 

Time dissemination experiments with ATS-1 

1968 

Initial STDN synchronization by GEOS-11 

May 1 968 | 

Initial NBS use of TV time comparisons between 
master clocks at Fort Collins and Boulder, 
Colorado. 

July 1970 

C-6and, two-way synchronization between NASA 
stations, through ATS-3. 
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6.1 FREQUENCY COMPARISON WITH VLF SIGNALS 


Radio stations NAA (17.8 kHz), NLK (18.6 kHz), WWVL (20 kHz), NSS 
(21.4 kHz) and WWVB (60 kHz) provide frequency reference coverage in 
the United States. Reception of a signal from any one of these sta- 
tions permits direct control of a local reference source or a cwnpari- 
son measurenent. The measurenent, or the control, refers the local 
standard frequency to the frequency standards at the National Bureau 
of Standards and at the Naval Observatory^^®*^^^. 

Receivers used for reception of these signals generally provide 

a measure of the phase difference between a local standard and a phase- 

locked or servo-controlled oscillator in the receiver. The phase- 

locked oscillator can be used as a frequency source or a long-term 

average of the phase difference between the two oscillators can be 

used to control or adjust the local standard. Propagation effects may 

q 

cause fluctuations of a few parts in 10 in the receiver oscillator 
frequency, but the 24 hour accuracy of a compa»^ison measurement is 
reported to be about 1 part in 10^^. 

The phase difference between a local standard and the received 
signal is recorded on a strip chart recorder and t>e change in phase 
difference over a 24 hour period is noted to deter line the frequency 
difference between the local and remote standards. As an example, 
consider the use of a signal from NAA at 17.8 kHz. Full scale deflec- 
tion on the chart is set for one cycle, 360 degrees, which corresponds 
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to 56.2 microseconds at the signal frequency. If then, the phase 
difference changes by 0.08 cycles over a 24 hour period, the corres- 
ponding time is (0.08)(56.2) * 4.496 microseconds. The fractional 
frequency offset is this phase difference, expressed in microseconds, 
divided by the period over which the measurement was made, in 
microseconds. 


Af , At ^ 4.496 

T IaTT io®(8.64 X 10^) 


5.2 X 10'^^ 


An increasing or decreasing phase indicates the local frequency is 
high or low, respectively. 


Propagati n time changes from day to night and is significantly 
variable at night. To reduce the effect of these changes in propaga- 
tion time, niqhtime measurements should be avoided. The 24 hour 
measurement of the example is not adversely affected by night effects 
if the readings are taken at the same time during dav hours. 

For good accuracy, initial and final phase readings are made when the 
entire propagation path is in the daylight. 


6.2 APPLICATIONS OF LORAN-C 


Much use has been made of Loran-C (100 kHz) signals for both time 
and frequency dissemination where high accuracy is required and where 
the ground wave signal strength is adequate^^®^ . In gep°>^al , an adequate 
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ground wave signal strength will be available within 1000 nautical 

miles of a Loran-C station if propagation is over land. This distance 

is increased to around 1400 nautical miles when the major length of 

the path is over sea water. Within this groundwave range, frequency 

13 

resolution between 5 and 10 parts in 10 has been reported. Third 
order servo control of a crystal standard that has high intrinsic sta- 
bility has been demonstrated to better than one part in 10 . 

Skywave measurements with Loran-C receiving equipment have proven 
feasible at distances of approximately 5000 miles^^^^. For this radius 
of coverage, all major land masses of the earth are covered by usable 
signals, with the exception of the continent of Antarctica^^^^. Timing 
accuracies for skywave measurements, are expected to be within 100 
microseconds even though the ionospheric conditions are not well known. 
With a one-time on-site calibration, to provide an estimate of the 
propagation delay, skywave timing synchronization is improved to about 
10 microseconds. 

In a typical remote instrumentation site, Loran-C might be used 
as suggested in Figure 13, for time and frequency synchronization. Two 
transmission paths are shown, one from the master station and another 
from the slave. The master station would be selected, if the signal 
strength is adequate, because of the more direct link to the USNO. At 
some locations, however, relative signal strengths may result in better 
measurements when the slave station signal is used. 
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STEERING INFORMATION 


k 
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Figure 13. - Loran C frequency and time services for an instrumentation site. 













Two local frequency 'Standards are shown at the remote site although 
only one is needed. An option of one or the other is suggested, or the 
use of both for redundancy. The exact configuration will vary, of 
course, with the situation at each site. WMV timing is used to set 
the local clock to an accuracy somewhat better than a second, perhaps 
10 milliseconds. A phase comparison is marie between one of the stan- 
dards and the time-averaged reference frequency from the Loran-C re- 
ceiver. The comparison could be visual, with a manual adjustment to 
the frequency standard. Alternatively, the phase comparator output 
can be used to control the frequency standard, establishing and main- 
taining coherence with the USNO standard. The f-^equency standard out- 
put is then distributed to the time code generator for precise time 
interval caliLrations. Time codes are recorded for time-correlation 
of test measurements. 

6.3 SYNCHRONIZATION THROUGH TELEVISION NETWORKS 

Frequency and time standards at the Newark Air Force Station (NAFS), 

United States Naval Observatory (USNO), and the National Bureau of Stan- 

(251 

dards have been compared, for a period of several years' , through 
signals from the nation-wide television networks. The comparison is 
made with relatively simple and inexpensive equipment which extracts 
the horizontal sync pulse for the tenth line of the odd field. These 
pulses are generated in New York City by reference to a rubidium fre- 
quency standard and are, therefore, relatively stable. 
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The composite TV signals are transmitted to the three stations over 
many microwave relays and are finally broadcast f.'om local VHP tele- 
vision stations. Operations at each of the three participating stations 
must be accurate to within 30 milliseconds, before the time comparison 
is made, to avoid the ambiguity that results from the repetitive nature 
of the television frame structure. At a preset time, to within one 
second, counters are started at each of the three statiors. The counters 
are then stopped at each station by the arrival <nd extraction of the 
odd field, tenth line sync pulse. Path delay calibration for each 
station and the number of high frequency clock cycles c.^nted at each 
station permits the calculation of differences in the time of occurrences 
of the one pulse per second clock pulses which started the counters. 

Measured time differences were within 6 microseconds over a period 
of more than a year. Typically, the mean standard deviation about the 
mean difference was about 0.5 microsecond over periods of one month. 

An improvement of about one order of magnitude is predicted for this 
technique through cooperative measures that could easily be taken at 
the originating studios, along the relay paths, and at the observers 
station. 
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6,4 SATELLITE DISSEMINATION OF FREQUENCY 
AND TIME INFORMATION 


Relatively simple equipment is described in Reference 26 for 
receiving time and frequency signals from the ATS-3 satellite. This 
method ot obtaining a long-term stability and accuracy .eference was 
considered. A survey of related literature revealed extensive acti- 
vity in this area. Two modes of operation have been practiced and 
experiments are continuing with both modes. One of these is passive 
and requires receiving equipment at the ground location, but no trans- 
mitting equipment. The other mode uses two way transmissions to ob- 
tain good measures of propagation delay. 

Synchronization by one-way links from the Tactical Communication 
Satellite (TACSAT) has been successfully completed. Five widely 
separated receiving stations in North and South America were time 
synchronized to within 150 microseconds. Analysis of the data revealed 
the possibility of improving the synchronization to 75 microseconds by 
a priori calibration of one s+ation with the master clock. Similar 
experiments with the Lincoln Experimental Satellite (LES-6) indicated 
an expected accuracy of 2b miv.'^oseconds, through the use ot these same 
methods . 

(271 

lime synchronization of a NASA network' , using a clock on-board 
GEOS-11 has also been demonstrated. In this operation, each station 
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recorded the difference between a locally generat'"-* ''lock pulse and a 
pulse received f>'*om the 6E0S-11 crystal -contr)l led clock. One of the 
stations was designated as the master and all other st‘'tions were 
synchronized to it by a conip=*r1son of the measured differences. In 
this way, the network clocks were synchronized, on an average, to 
within 25 microseconds of the master clock at the c ^ignated master 
station. 

Two-way satellite synchronization experiments have demonstrated 
a greater accuracy than any other satellite technique, competing with 
the most accurate methods now under investigation. A portable cesium 
clock was used as a basis for evaluation in suc'i a tv station experi- 
ment. Ihe time difference of the ground station clocks ; as measured 
over the two-way satellite links, differed from that measured with a 
portable cesium clock by less than 50 na .oseconds. The sii;.tf;cant 
improvement over one-way link measurements is in the propagation time 
calibration that is in the .need for transmitting equipment at each 
station to be calibrated. When the relative positions of network 
stations are well known, however, and when these stations have some- 
what similar ray paths, excellent synchronization may be oossible with 
transmitting equipment at only one or two stations. 

Satellite dissemination methods, when compared with VLF and LF 
dissemination, appear to require more expensive equipment ano .-re rot 
as well developed. In addition, the techniques app'‘ar to be orirnlod 
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toward an onphasis on time dissemination. This is somewhat of a dis- 
advantage, even though the time corrections ultimately correct the 
frequency standard which drives the clock. 
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7.0 AN ERROR MODEL FOR FREQUENCY STANDARD STABILITY 


A mode's of frequency standard stability is needed for two pur- 
poses. First, a model is needed for use in the evaluation of Doppler 
extractor performance and the Doppler measurement configuration that 
is to be used for the Shuttle. Then, a second model is needed, simp- 
lified to reduce da*:a processing requirements, for use in the naviga- 
tion program on-board the Shuttle. This simplified model might in- 
clude any effects of the ground station stability that are large 
enough to introduce significant navigation errors. 

The model should include the elements on the r’ght of Equation 2 
and additional elements to account for environmental effects. System- 
atic frequency deviations snould be included as functions of their 
causative variables. The noise term w’ll be increased by some of 
the environmental effects such as vibration. To include these effects, 
the variance value which represents each of them will be summed with 
the square of n^. 

f(t) = . V . * ["/ ^ 

Systematic effects of temperature, line voltage and acceleration 
should be represented by terms in the model. Residuals from these 
terms should be included as random effects, along with the random 
effect of vibration. C„ray noise, pi eke i up by radiation and conduction. 
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may also be significantly large. If so, it would appear as a random 
effect. It should be noted that the random terms may have amplitudes 
that vary with flight operations. Vibration is an obvious example. 
Stray noise may be another. 


7-2 



8.0 CONCLUSIONS AND RECOMMENDATIONS 


Frequency stability literature covers both frequency and time 
domain descriptions. Because of the range of time intervals over w*ich 
Doppler counts are planned in the Shuttle application, time domain 
descriptions and measur^ents seem preferable for evaluation of the 
Shuttle Doppler measurements and the equipment used in making them. 

It may be desirable, however, to make a relatively small number of 
frequency domain measurements. 

Methods already in use, using period counting techniques in the 
time domain, seem satisfactory and are similar to methods that have 
been reported in the literature. These methods are well suited for 
short-term stability measurements and can be used over a period of 
time to collect data from which long-term drifts can oe calculated. 

The drifts are not calibrated, however, on the frequency axis. Such 
calibration appears feasible through the use of a VLF/LF receiver. 

By this means, the calibration can be referred to frequency stan- 
dards at the National Bureau of Standards and the Naval Observatory. 

The procedure used would also provide a simpler, mure direct and more 
accurate measure of long-term stability. 

A survey of literature and commercial specifications shows the 
performance of a ruggedized crystal oscllator to be at least as good, 
in almost all respects as a good commercial crystal standard. In 
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fact, performance specifications for a crystal standard that is 
designed for satellite service are generally better than those for 
comnercial standards. This same trend is not so noticeable for rubi 
dium standards, for which some degradation appears to result from an 
attempt to reduce size and weight. 

Available information on environmental effects appears to be 
adequate for a first model approximation. Tests on the Shuttle 
oscillator will be needed, however, to refine the model. A possible 
exception to this conclusion regarding the adequacy of available in- 
formation exists in the area of radiation effects. A radiation 
specification for the Shuttle is not available at this time. 
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APPENDIX A 

DATA FOR CRYSTAL OSCILLATOR 
FREQUENCY STANDARDS 



The selection of transmitting and receiving frequency source types 
for making one-way Doppler measurements on-board the Shuttle Orbiter 
have been made. A significant factor in each selection was the stability 
characteristics of the various contenders. Simplicity was also a signi- 
ficant factor in the selection of a crystal oscillator for use on-board 
the Shuttle because of the decreased cost, size, and weight whic''. accom- 
pany designs of lower complexity. As a general rule, simplicity also 
enhances reliability. 

Circuit details of the Shuttle frequency reference oscillator are 
not yet known. In some respects, however, it is expected to be like 
the exemplary circuit of Figure A-1. This illustration includes most 
of the elements which are generally found in crystal oscillators that 
are intended for use as stable frequency sources. 

The double oven maintains a constant crystal temperature over a 
specified operating range, to a tolerance of about 0.01°C. Heater 
current is not cycled off and on abruptly, but is smoothly controlled. 
Heating current is proportional to the difference between the actual 
and desired temperatures. Although reasonably good supply voltage 
regulation is expected, addi..iona1 regulation is usually provided 
for all active elements. Automatic gai;- control (AGC) is used to 
reduce the sensitivity of the oscillator frequency to signal level 
changes. Both the crystal and the active elements undergo parameter 
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Figure A-1. — A functional illustration of an exemplary crystal oscillator. 




changes as signal levels vary. Another crystal Is sometimes used as a 
filter, following amplification which brings the signal level above the 
very low level of the oscillator. (Low crystal drive levels are used 
to improve stability.) Filtering with a second crystal reouces the 
noise sideband levels, particularly at frequency offsets greater than 
10 Hz. A buffer amplifier isolates the oscillator from changes in 
loading, reducing the sensitivity of the oscillator frequency to load 
resistance. 

All of the elements in Figure A-1 are not always included in a 
good oscillator design. The elements that are included in a particular 
design will depend, to some extent, on tiie application. For example, 
crystal filtering would not be necessary in applications where noise 
modulation at frequencies above 10 Hz is relatively unimportant. 

Characteristics for oscillators of relatively good design are 
sunmiarized in Table A- I. Only the FE 1800D is specified for space- 
borne applications. For this reason, information on performance in a 
space environment is not generally available. The FE 1800D, however, 
is specified to perform as well or better than other units in almost 
every case. It might be concluded, then, ti«at the ruqgedized design 
does not degrade performance. Oscillator cha^actei sties for several 
test instruments are listed in table A-II. 
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TABLE A-I - COMMERCIAL CRYSTAL S-ANDARD CHARACTERISTICS 
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TABLE A-II. - TEST EQUIPMENT OSCILLATOR CHARACTERISTICS 


















APPENDIX B 

ATOMIC FREQUENCY STANDARDS 


/ 



1.0 APPLICATIONS OF ATOMIC RESONANCE PHENOMENA 


Two basic techniques have been used, in the development of fre- 
quency standards, which rely on the intrinsic precision of atomic 
phenomena. These two techniques are illustrated in Figure B-1 and, 
as noted there, can be distinguished by the descriptive terminology, 
“passive resonator", and "active oscillator". In the passive techni- 

an atomic device is used as a high-Q resonant circuit component 
to control the frequency of a crystal oscillator. Rubidium gas cells, 
cesium beam tubes, and thallium beam tubes have been used as the reso- 
nant elements in frequency standards of this type. In the active tech- 
nique, a crystal oscillator is phase locked to a self-excited, atomic 
oscillator. Hydrogen makers, ammonia masers, and optically pumped 
rubidium masers have been used as atomic oscillators in this type of 
frequency standard. 

Atomic resonance characteristics are centered above 1000 MHz for 
both passive and active devices which have been used in the design of 
commercial and laboratory frequency standards. The crystal oscilla- 
tors that are used in these standards generally operate at 5 or 10 
MHz, but sometimes as high as 100 MHz. Passive rubidium cells, passive 
cesium beam tubes > and hydrogen maser oscillators have proven most 
useful in commercial and laboratory frequency standards of recent 
design. 


preceding ?AGE RLANK KG'? FILMED 



ATOMIC 

RESONATOR 




ATOMIC 

OSCILLATOR 


Figure B-1. 



(a) Passive atomic resv.nator 



(b) Active atomic oscillator 


— Essential elements for two atomic frequency 
standard designs. 
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AtcHnic frequency standard characteristics are ccunpared in Table 
6-1. Characteristics are also listed there for a good crystal oscilla- 
tor. In fact, the values listed in Table B-I are generally among the 
best that can be found. The minimum flicker fluctuation is an indi- 
cation of the best stability that might be achieved with a given type 
of standard and is independent of the averaging period. 

The long-term stability of atomic standards is an important 
advantage, compared to a crystal standard. The higher cost of atomic 
standards is a disadvantage, of course. For spacebome applications, 
power and weight favor the crystal standard. The rubidium performance 
figures in Table B-I are for one of the larger units, but rubidium 
standards are seen to compete with the crystal standard in size. 
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TABLE B-I - A COMPARISON OF FREQUENCY STANDARD? 
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2.0 RUBIDIUM FREQUENCY STANDARDS 


Rubidium gas cells have been used in both active and passive 
frequency standard designs. As an active frequency standard element, 
the rubidium maser is relatively simple and compact, when compared to 
the hydrogen and ammonia standards. Light intensity and frequency 
pulling effects introduce drift and result in flicker noise levels 
that are somewhat higher than other masers. 

Commercial designs that use rubidium gas cells have favored the 
passive method of frequency control. A summary of commerical rubi- 
dium standard characteristics is presented in Table B-II. Data for 
Table B-II have been taken from specifications provided by the four 
companies whose products are represented there. The FRK units have 

(33 34 35) 

been designed and tested for space applications. Test results' * ’ ' 

show the ultimate long-term drift to be better than one part in 10^^ 

per month after a period of stabilization. In fact, the drift was 

found to be "essentially zero" after an initial drift rate of about 
12 

2.5 parts in 10 per day. 

Sensitivity to acceleration at g-forces just above zero was found 

-12 -12 
to be about 8 x 10 per g. The sensitivity dropped to 4 x 10 /g 

at Ig, however, and remained at this value through the maximum level 
(lOg) of the test. This sensitivity would normally imply good per- 
formance under vibration, if mechanical resonances were adequately 
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TABLE B-II - CHARACTERISTICS OF RUBIDIUM FREQUENCY STANDARDS 
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(See next page) 














TABLE B-II (Continued) 


NOTES: 

(1) lx 10”^^ for any orientation in the earth's fie'd. 

(2) Also, curve 1 of MIL-E-5400, with isolators. 

(3) Also, MIL*E-5400 (30 g). 

(4) 0.005 X 10 for 1 gauss change. 

-Q 

(5) 0.05 X 10 for altitudes from sea level to 40,000 feet. 


B-8 



suppressed. Performance under vibration was not acceptable, however, 
because of equipment failures. Corrective designs were not believed 
to be extensive, but no further test data are available. 

Temperature sensitivity varied for two units which were subjected 
to a range of temperatures from 5 to 40°C. For one unit, the variation 
was essentially linear through a total variation of about 1 part in 
10^®. Over a selected 5° range, this unit had a 0.4 x 10”^^/°C devia- 
tion. Performance o the second unit was more erratic, with a high 
variance around a somewhat parabolic characteristic. Over a selected 
5° range, it had a 0.3 x 10"^®/°C deviation. 

The sensitivity to ambient pressure was measured between standard 

-6 

atmospheric conditions and a pressure of 10 Torr. (A Torr is equiva- 

3 

lent to 1.333 x 10 dynes/per square centimeter or 1 millimeter of mer- 
cury. Standard almiospheric pressure is 760 millimeters of mercury. 

The international unit of pressure, the pascal, is equivalent to 10 
dynes per square centimeter or one newton per square meter.) Over this 
range, the sensitivity was negative. The frequency increased 12.3 parts 
in 10^® for the decrease in pressure. 

Sensitivity to radiation was relatively low. It was reasoned that 
this low sensitivity resulted from counter-acting effects. An expected 
change in the crystal was offset by a change in control voltage at the 
time of irradiation. A 300 rad dose was applied over a period of 10 
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hours, at a constant rate of 0.5 rads/minute. (A rad Is equivalent 
to 10 joules per kilogram.) A change of -1.5 parts In 10 was 
observed when the Irradiation was begun. No further change was 
noted, even for a second Irradiation. The effect remained evident 
for a period of 5 days. 

Fractional frequency fluctuation over a period of one second was 
measured to be about 2 parts In 10^^. This decreased to about 1 part 

IP 

In 10 for a 100 second period. 

A comparison of the rubidium standard stabilization or warmup 
time with that of a crystal standard reveals what might be a slgnlfl^ 
cant advantage. One of the rubidium standards noted In Table B-II 
Is specified to be within 2 parts In 10^^ of Its operating frequency 
after a 10 minute warmup. Crystal standard warmup time Is nominally 
one order of magnitude greater thai! this, at best. 
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3.0 CESIUM-BEAM STANDARDS 


Commercial cesium beam frequency and time standards use a cesii"r> 
tube which has an electrical resonance characteristic. The phase 
response of this resonant element shifts the phase of a control signal 
that is used in a servo loop to control the frequency of a crystal 
oscillator. 

As noted in Figure B-2, the output frequency from a crystal oscilla- 
tor is passed through a synthesizer and a multiplier chain to a mixer. 

A phase modulated sample from the crystal oscillator is also applied 
to the mixer, producing a phase modulated signal at, or near, the 
resonant frequency of the cesium beam tube. If this signal differs 
in frequency from the resonant frequency of the cesium beam tube, the 
beam current of the tube will contain a component at the frequency of 
the control signal. The phase of this componen;: of current depends 
on the direction of the difference between the cesium tube resonant 
frequency and the synthesized frequency from the crystal oscillator. 

The component of cesium tube current at the control frequency is 
selected by a tuned amplifier and passed to a phase detector. A sample 
of the control frequency is also aoplied to the phase detector, directly 
from the control oscillator, as a reference for comparison with the 
beam tube current component. The phase difference between these two 
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Figure B-2. - A typical cesium beam frequency standard. 






signals produces an error signal, e (</>), at the output of the phase 
detector. An integrator is used to smooth the error signal and this 
smoothed signal, J*e(<#>)dt, corrects the crystal oscillator frequency. 

The extremely stable operation of the cesium beam standard is due 
to the very high Q of the resonant response of the tube. A half- 
amplitude bandwidth of 550 Hz has been measured, indicating a Q of 
18 million, at the center frequency of over 9000 megahertz. 

Cesium beam time standards have been transported by airplane and 
automobile to points throughout the world. This endeavor has demon- 
strated the rugged and dependable qualities rf a "flying clock" and 
the practicality this method of comparing widely separated time stan- 
dards. The reliability of laboratory models of the cesium beam stan- 
dard is given as better than a 15,000 hour MTBF. The cesium beam tube 
is specified to i«*ve a 30,000 hour MTBF. 


N*SA-JSC 
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